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1 Introduction 

Recall that a Schrodinger operator on a n-dimensional Riemannian manifold (M, g) is a second 
order linear differential operator of the form 

h = ~a + u, 

where A is the Laplace-Beltrami operator and U is the potential function for the physical system 
under consideration. A question of fundamental interest in quantum mechanics is to construct 
eigenfunctions ip of the Schrodinger operator TCq. One approach to this problem, which is based 
on the representation theory of Lie algebras, is to consider Schrodinger operators TCq which are 
quasi-exactly solvable, in a sense that will be defined below. 

We begin by considering the case of a general linear second-order differential operator 7i, 
given in local coordinates by 

n n 

n=J2 Aijd i d i + B ' di + a 

i,j=l i=l 

The operator TL is said to be Lie algebraic if it is an element of the universal enveloping algebra 
of q, a finite dimensional Lie algebra of first order differential operators. More explicitly, 

m m 

n=J2 C -bT a T b + c a r a + C , (1) 

a,b=l a=l 

where 



T a = v a + rj a , l<a<m, (2) 

is a basis or rj. In the operators v m are vector fields, and rji, . . . , rj m are multiplication 

operators. A Lie algebra g of first order differential operators is said to be quasi-exactly solvable 
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if one can find explicitly a finite dimensional g-module N of smooth functions, i.e., if J\f = 
{h , . . . ,h r } with T a (N) C M for all 1 < a < m. A Lie algebraic operator 7i, is said to be 
quasi- exactly solvable if it lies in the universal enveloping algebra of a quasi-exactly solvable 
Lie algebra of first order differential operators. Obviously, 7i{M) C J\f, i.e. the module M will 
be fixed by the operator 7i. Moreover, if the functions contained in the module M are square 
integrable with respect to the Riemannian measure ^Jgdx 1 ■ ■ ■ dx n , where g is the determinant of 
the covariant metric, the operator 7i is said to be a normalizable quasi-exactly solvable operator. 

We can see from the above definitions that the formal eigenvalue problem for quasi-exactly 
solvable Schrodinger operators can be solved partially by elementary linear algebraic methods. 
Indeed, the operator H is self-adjoint with respect to the inner product associated to the standard 
measure, therefore the restriction of 7i to the finite dimensional module M is a Hermitian finite 
dimensional linear operator. Thus, one can in principle, compute r = dim(AA) eigenvalues of 7i, 
counting multiplicities, by diagonalizing the r x r matrix representing TC in a basis of M. 

It seems that the concept of a "spectrum generating algebra" was first introduced by Goshen 
and Lipkin in [15] in 1959. However, this paper did not seem to have been be noticed by the 
community and, ten years later, spectrum generating algebras were independently rediscovered 
by two groups of physicists, see [2] and [5]. Their work was an impetus for further research in 
this area as one can see by browsing in the two volume set of reprints [3] and the conference 
proceedings [16] . A survey of the history and the contribution papers related to the spectrum 
generating algebras is given in the review paper of Bohm and Ne'eman, which appears at the 
beginning of [1]. In the early 1980's, Iachello, Levine, Alhassid, Giirsey and collaborators ex- 
hibited applications of spectrum generating algebras to molecular spectroscopy; a survey of the 
theory and applications is given in the book [T7]. In these applications, both nuclear and spec- 
troscopic, the relevant Hamiltonian is a Lie algebraic operator in the sense described previously. 
Finally, the analysis of a new class of Schrodinger operators, the quasi-exactly solvable class, 
was initiated in late 1980's by Shifman, Turbiner and Ushveridze, see |23 |. 124 1 125] . A survey of 
the theory and applications of quasi-exactly solvable systems in physics is given in [26] . 

There exists a complete classification of quasi-exactly solvable Schrodinger operators in one 
dimension. In two dimensions, this classification is in principle complete. Indeed, all the Lie al- 
gebraic linear differential operators for which the formal spectral problem is solvable are known. 
The question of determining if the operator is equivalent to a Schrodinger operator will be 
discussed below. The main contribution of our paper is to extend these results to three dimen- 
sions by giving a partial classification of the quasi-exactly solvable Lie algebras of first order 
differential operator in three variables, and showing how this can be applied to the construction 
of new quasi-exactly solvable Schrodinger operators in three dimensions. Our work is based on 
the classification of finite dimensional Lie algebras of vector fields in three dimensions begun 
by Lie in [18] and completed by Amaldi in pQ. 

Let jW) be the contravariant metric of the manifold M in a local coordinate chart, g its 
determinant and g the determinant of the covariant metric. In that setting, a Schrodinger 
operator reads locally as 



1 n 

«o = -i£ 
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^d tJ + d, l (g^)d J - 9 -^d J 



+ U. 



Recall that a quasi-exactly solvable second order operator is not, in general, a Schrodinger 
operator. However this operator might be equivalent to a Schrodinger operator in a way that 
preserves the formal spectral properties of the operators under consideration. The appropriate 
notion of equivalence, which will be used throughout our work, is the following. Two differential 
operators are locally equivalent if there is a gauge transformation TL — ► fiTCfj, -1 , with gauge 
factor (j, = e x , and a change of variables relating one to the other. In principle, it is possible to 
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verify if a general second order differential operator 7i is equivalent to a Schrodinger operator 
with respect to this notion of equivalence. Indeed, every second order linear differential operator 
can be given locally by 

.. n n 

H = -2 T.^dij + J^h'di + U. 

i,j=l i 

If the contravariant tensor g^> is non-degenerate, that is if g does not vanish, the operator can 
be expressed as 

H = -^A + V + U, (3) 

where V = b l di is a vector field. For this operator to be locally equivalent to a Schrodinger 
operator, the vector field V has to be a gradient vector field with respect to the metric 
Locally, this will be the case if and only if uj = gijVdx 1 , the one form associated to V, is closed. 
For this reason, this condition is named the closure condition. Note that if V = V(A), the gauge 
factor is given by e~5. 

Given an operator of the form ([I]), the closure conditions can be easily verified provided the 
contravariant metric g^ is non-degenerate. Indeed, these conditions can be written as algebraic 
constraints on the coefficients C a b and C c , and are the Frobenius compatibility conditions for 
an overdetermined system that will be described later. 

An important point to keep in mind is that the class of quasi-exactly solvable operators is 
invariant under local equivalence. Indeed, suppose Ti is a quasi-exactly solvable operator which 
is gauge equivalent to an other operator 7i$ under the rescaling [i. If TC lies in the universal 
enveloping algebra of g, whose g-module is Af, one can easily show that TLq is quasi-exactly 
solvable with respect to the finite dimensional Lie algebra 

Q = f i-Q-^ 1 = {fi-T-pT 1 | Teg} 

which is isomorphic to g and posses the finite-dimensional g-module 

N = n-j\[ = {n-h\he M}. 

Note however that the gauge factor is not necessarily unitary. Thus, a gauge transformation does 
not necessarily preserve the normalizability property of the functions in Af. Therefore, the class 
of normalizable quasi-exactly solvable operators is not invariant under our local equivalence. 
We now give an example of a normalizable quasi-exactly solvable Schrodinger operator in three 
variables. 

Example 1. In this example we consider the quasi-exactly solvable Lie algebra g = s[(2) xst(2) x 
5l{z). With the standard notation p = q = and r = J-j, this Lie algebra representation 
can be spanned by the following first order differential operators 

T 1 = p, T 2 = xp, T 3 = x 2 p — x, T 4 = q, T 5 = yq, 
T 6 = y 2 q -y, T 7 = r, T 8 = zr, T 9 = z 2 r - z. 

Then, the finite dimensional module of smooth functions 

Nm x m y m z ■= {x % y j z k \ < i < m x , < j < m y , < k < m z } 

is a g-module provided m x = m y = m z = 1. With the following choice of coefficients, one 
constructs the quasi-exactly solvable operator 

-2H = (T 1 ) 2 + (T 2 ) 2 + 2[(T 3 ) 2 + (T 4 ) 2 + 2(T 5 ) 2 + 2(T 6 ) 2 ] + (T 7 ) 2 + (T 8 ) 2 
+ 2(T 9 ) 2 + {T 1 , T 3 } + {T 7 , T 9 } - 2T 4 - 2T 5 - 4T 6 - 8, 
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where {T a , T b } = T a (T b ) +T b (T a ) . The induced contravariant metric associated to this operator 
is computed to be the following positive definite matrix 

/ (x 2 + l) 2 {x 2 + l){z 2 + l) \ 

g m =\ 2/ 4 + 4y 2 + l , (4) 

V (x 2 + l)(z 2 + 1) 2(z 2 + l) 2 / 

whose determinant is g = (x 2 + l) 2 (y 4 + 4y 2 + l)(z 2 + l) 2 . Then, with respect to this non- 
degenerate metric, the operator TL can also be described as 

-2H = A + V + U, 

where V = — 2(x 3 + x + z + zx 2 )p — 2(2y + y 3 )q — 2(2z 3 + 2z + x + xz 2 )r. It is not hard to verify, 
always with respect to the metric (jH , that the first order term V is the gradient of the function 
A = — ln(x 2 + 1) — l/21n(y 4 + 4y 2 + 1) — ln(z 2 + 1). Hence, by considering the gauge factor 



A 



/2 



M = e2 = ( X 2 + 1)-V2 (y 4 + Ay 2 + 1} -l/4 (z 2 + 

the operator H is gauge equivalent to a Schrodinger operator 
-2W = A + U, 

were the potential is a rational function of y. Furthermore, it is not hard to show that, after the 
gauge transformation, the functions in J\fm = { [i-x t y i ' z k \ < i, j, k < 1} are square integrable 
with respect to yjgdxdydz. Recall here that g is the determinant of the covariant metric, hence 
■yjg = fj, 2 . Thus, for i, j, k either or 1, one can use Fubini's theorem to decompose the integral 

r w;*l 2 „ 2 ^,,^ = fff x 2l y 2 ^z 2k 



(fix l y J z ) ix dxdydz = , 2 - , 2 , 4 - 2 - v 2 2 dxdydz, 

i JJJrs (a: 2 + l) 2 (y 4 + 4y^ + 1)(^ + l) 2 

into the product of three finite integrals in one variable. Consequently the operator T~Lq is 
a normalizable quasi-exactly solvable Schrodinger operator and it is possible to compute eight 
eigenfunctions by diagonalizing the matrix obtained by restricting Ti to M. For this operator, 
one gets two eigenvalues, —3 and 1, both of multiplicity four. The eight eigenfunctions associated 
to these two eigenvalues are respectively, 

"0-3,1 = -1 + xz, 0-3,2 = y ~ xyz, V-3,3 = xy + yz, il>- 3 ^=x + z, 
■0i,i =y + xyz, ip lt2 = -x + z, Vi.3 = -xy + yz, ip 1A = l + xz. 

Finally, one gets eight eigenfunctions of the Schrodinger operator 7Yo by scaling each of these 
functions by the gauge factor //. 

In general, there is no a-priori method for testing whether a given differential operator is Lie 
algebraic or quasi-exactly solvable. However, one can try to perform a classification of these 
operators under local equivalence using the four-step general method of classification described 
by Gonzalez-Lopez, Kamran and Olver in [llj . 

The first step toward the classification of normalizable quasi-exactly solvable Schrodinger 
operators is to classify the finite dimensional Lie algebras of first order differential operators 
up to diffeomorphism and rescaling. Then, the task is to determine which of these equivalence 
classes admit a finite dimensional g-module N of smooth functions. Then, from the quasi- 
exactly solvable Lie algebras found in the second step, one can construct second order differential 
operator as described in ([1]) from any choice of coefficients C a b, C c , and Co- The third step 
consists to determine which of these operators are equivalent to Schrodinger operators and this 
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can be performed by verifying the closure condition. Finally, the last step in this classification 
problem is to check if the functions contained in the g-module N are square integrable. 

As mentioned previously, the entire classification has been established in one dimension. In 
the scope of the first two steps, every quasi-exactly solvable Lie algebra is locally equivalent to 
a subalgebra of the Lie algebra 

c i d d 2 d A 

Q n = Span <^ — , x— , x nz , 1 } , 

[ ox ox ox J 

where re is a non negative integer, see [13] for more details. Then, once a second order differential 
operator is constructed, since all one forms are closed in one dimension, such operator will 
always be equivalent to a Schrodinger operator, reducing the third step to a trivial step. Finally 
Gonzalez-Lopez, Kamran and Olver determined in [TO] necessary and sufficient conditions for 
the normalizability of the eigenfunctions of the quasi-exacly solvable Schrodinger operators. 

In two dimensions, the first two steps of the classification problem were determined by the 
same authors in [T2J and [14] . Based upon Lie's classification of Lie algebras of vector fields, 
see [18] . a complete classification of the quasi-exactly solvable Lie algebras q of first order 
differential operators, together with their finite dimensional g-modules, was completed. The 
case of two complex variables is discussed in the first two papers while the third paper completed 
the classification by considering operators on two real variables. However, the last two steps 
are not yet completed but a wide variety of normalizable quasi-exactly solvable Schrodinger 
operators has been exhibited, see for instance fl~3] and [14] . 

In the next section, a partial classification of quasi-exactly solvable Lie algebras of first order 
differential operators in three dimensions is given. While these two first steps were successfully 
completed in one and two dimensions, only part of this work is now done in three dimensions. 
However, these new quasi-exactly solvable Lie algebras can be used to seek new quasi-exactly 
solvable Schrodinger operators in three dimensional space. The last section of this paper is de- 
voted to the description of new quasi-exactly solvable Schrodinger operators in three dimensions. 
Eigenvalues are also computed for two families of Schrodinger operators. These eigenvalues are 
part of the spectrum of the operators and their eigenfunctions, together with their nodal sur- 
faces, are exhibited. In addition, a connection is made between the separability theorem proved 
in [6] and the quasi-exactly solvable Schrodinger operators on flat manifold. The quasi-exactly 
solvable models obtained in our paper are new as far as we can tell. In particular they are 
not part of the list of multi-dimensional quasi-exactly solvable models obtained in [26] by the 
method of inverse separation of variables. 

2 Classification of quasi-exactly solvable Lie algebras 
of first order differential operators 

2.1 Lie algebras of first order differential operators 

Our goal in this section is to give a partial classification of quasi-exactly solvable Lie algebras of 
first order differential operators in three dimensions. A first step toward this goal is to obtain 
a classification of the finite dimensional Lie algebras of first order differential operators. 
After this is done, the next step is to impose the existence of an explicit finite dimensional g- 
module J\f of smooth functions. To this end, we will first summarize the basic theory underlying 
the classification of Lie algebras of first order differential operators. 

For M an re-dimensional manifold, we denote by ^"(M) the space of smooth real-valued 
functions and V(M) the Lie algebra of vector fields on M. The space ^-"(M) form a V(M)- 
module under the usual derivation rj — > v(rj), where v is a vector field in V(M) and r/ a function 
in JF(M). The Lie algebra of first order differential operators P 1 (M) can be described as 
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a semidirect product of these two spaces, 2? 1 (M) = V(M) x JF(M). Indeed, each element T 
in P 1 (M) can be written into a sum T = v + r\ and the Lie bracket is given by 

[T\T 2 ] = [v l ,v 2 ] +v 1 (rj 2 ) -v 2 {t] 1 ), where T = v i + rf G ^(M). (5) 

Note that the space ^"(M) is also a 'D 1 (M)-module with T(£) = v(C) +rj ■ Consequently, any 
finite dimensional Lie algebra of first order differential operators q can be written as 

T 1 =V 1 + rj l , T s = v s +r] s^ T S+1 =C\ _ ^ TS+r = c r^ (g) 

where u , . . . ,v s are linearly independent vector fields spanning fj C V(M), a s-dimensional Lie 
algebra and where the functions . . . , £ r act as multiplication operators and span .A/f C ^-"(M) 
a finite dimensional f)-module. Note that restrictions need to be imposed to the functions rf 
for g to be a Lie algebra. Indeed, without the cohomological conditions that will be described 
below, the Lie bracket given in ([5]) does not necessarily return an element in the Lie algebra q. 

For T = v + 7], we define a 1-cochain F : \) — > ^-"(M) by the linear map (F; v) = 7]. Since 
any function ( £ M can be added to T without changing the Lie algebra q, this map is not 
well defined. To deal with this issue, we should therefore interpret F as a JF(M)/-M-vahied 
1-cochain. Thus, from the Lie bracket given in (0), it is straightforward to see that g is a Lie 
algebra if and only if the 1-cochain F satisfies the bilinear identity 

v l (F;v j ) -v^F-v 1 ) - (F; [v\v j ]) G M, v\v j et). (7) 

In terms of Lie algebra cohomology, this condition can be restated as follow, {5iF; v l , v 3 ) £ M. 
for all v l , n> in f), i.e. F is a •7 r (M)/.M-valued 1-cocycle on f). (See [7] for a detailed description 
of Lie algebra cohomology.) 

This classification of Lie algebras of first order differential operators would not be complete 
without considering the local equivalences between the Lie algebras. Indeed, if a gauge trans- 
formation with gauge factor = e\ is performed on an operator T = v + r\ in q, the resulting 
differential operator T = e x ■ T ■ e~ A = v + r/ — v(X) will only differ from T by the addition of 
a multiplication operator v(X). Again, this can be expressed in cohomological terms. Indeed, 
under the O-coboundary map So : f) — > ^"(M)/^ defined by (5o\;v) = v(X), the multiplication 
factor v(X) can be interpreted as the image, or the O-coboundary, of the function A. Hence, 
combining these two observations, it is possible to conclude that the map F is an element in 
H 1 ^, ^(M.) / A4) = ker 5i/Im<5(). Thus, if two differential operators q and q are equivalent with 
respect to a change of variables ip and a gauge transformation given by fi = e A , these two 
operators will correspond to equivalent triples (f), .A/f , [F]), and (f), Ai, [F]), where fj = (/?*(()), 
M = tp*(A4), and F = (p* o F o cp^ 1 + 5qX. This is summarized in the following theorem. 

Theorem 1. There is a one to one correspondence between equivalence classes of finite dimen- 
sional Lie algebras q of first order differential operators on M and equivalence classes of triples 
(t),M,[F\), where 

1) f) is a finite dimensional Lie algebra of vector fields; 

2) M is a finite dimensional ^-module of functions; 

3) [F] is a cohomology class in H 1 (l),F(M)/M). 

Hence the general classification of finite dimensional Lie algebras of first order differential 
operators q can be bring down to the classification of triples {\),M, [F]) under local changes of 
variables. 

In three dimensions, a complete local classification of the finite dimensional Lie algebras 
of vector fields f) has been established by Lie in [18] and Amaldi in p]. Lie's classification 
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distinguishes between the imprimitive Lie algebras, for which their exists an invariant foliation 
of the manifold, and the primitive Lie algebras, for which no such foliation exists. Lie's work 
gives a description of the eight different classes of primitive Lie algebras and, based under the 
possible foliations of the manifold, the imprimitive Lie algebras are subdivided into the following 
three types: 

I. The manifold admits locally an invariant foliation by surfaces that does not decompose 
into a foliation by curves. 

II. The manifold admits locally an invariant foliation by curves not contained in a foliation 
by surfaces. 

III. The manifold admits locally an invariant foliation by surfaces that does decompose into 
a foliation by curves. 

Observe that these three types are not necessarily exclusive. For instance, the Lie algebra 
f) = {p, q, xq, xp — yq, yp, r} belongs to the first two types. The underlying manifold M 3 admits 
a first indecomposable foliation by planes A := {z = constant} and also admits an second 
invariant foliation by straight lines $ := {x = constant} n {y = constant} not contained in any 
invariant surfaces. Lie classified the algebras of type I and II, giving respectively twelve and 
twenty-one different classes of Lie algebras. Few years latter, the 103 classes of Lie algebras of 
the third type were exhibited by Amaldi. 

The number of finite dimensional Lie algebras of vector fields \] is large and it did not seem 
reasonable to consider all the 154 classes. For this first classification attempt, we have chosen to 
focus on the algebras which seem promising in our aim to construct new quasi-exactly solvable 
Schrodinger operators. The selection was made upon the following criteria. 

We first narrowed our choice based on the results given in [6]; provided the Lie algebra q 
is imprimitive and its invariant foliation consists of surfaces, one can show, adding some other 
hypothesis on the metric induced, that a Lie algebraic Schrodinger operator generated by q 
separates partially in either Cartesian, cylindrical or spherical coordinates. Since such algebras 
are good candidates for generating interesting quasi-exactly solvable Schrodinger operators, we 
restricted our search on the type I and type III imprimitive algebras. In this paper, the clas- 
sification of the twelve type I Lie algebras is entirely performed while, for the type III Lie 
algebras, we focused on some of the most general Lie algebras. Since the induced metric g^' 
needs to be non-degenerate, the type III Lie algebras involving only one or two of the three 
partial derivatives were discarded. Finally we selected our algebras among those that contain 
other type III algebras as subalgebras. 

2.2 Classification of Lie algebras of first order differential operators 

Using the equivalence between the Lie algebras of first order differential operators and triples 
(f), Ai, [F]), it is possible to determine the Lie algebras q from the selected Lie algebras of vector 
fields f). But first, recall that the second step in the classification of quasi-exactly Schrodinger 
operators is to determine which of these Lie algebras of first order differential operators are 
quasi-exactly solvable. It is not hard to see that if is quasi-exactly solvable with non trivial 
fixed module M, the Lie algebra is finite dimensional if and only if M is the module of constant 
functions, see [12] for details. Therefore, instead of working on the general classification of Lie 
algebras of first order differential operators 0, we will restrict our work to the equivalence classes 
of triples (h, {1}, [F]). Thus, for each of the selected Lie algebras f), we first seek for the possible 
cohomology classes, [F] in H 1 (i},J r ('M.)/{l}). Once this is done, in the scope of the second step, 
it will be left to find if there exists an explicit finite dimensional 0-module M, where is the 
Lie algebra equivalent to the triple (h, {1}, [F]). Note that, as in lower dimensions, the existence 
of a nontrivial module M will impose a "quantization" condition on [F]. Indeed, for each of the 
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Lie algebras worked out in this paper, the possible values for the functions in [F] can only be 
taken in a discrete set. For detailed results related to the quantization of cohomology, see [8] 
and [22]. 

2.2.1 Classification of the cohomology classes [F] in JF(M)/{1}) 

To determine the possible cohomology classes, we first start with [F] as general as possible. For 
every v in the Lie algebra f), we denote the value of the 1-cocycle (F; v) by r] v , and r] v can be any 
function in .F(M)/{1}. Our aim is to find the most general 1-cocycle F, that is the most general 
functions rj v , satisfying the restrictions imposed by the 1-cocycle conditions ([7]). Then, using the 
O-coboundary map, we try to describe the class [F] with representatives rj v as simple as possible. 
Finally, if {v , . . . ,v r } is a basis for t), the set {v 1 + r) v i , . . . , v r + rj v r} will be a basis for the 
Lie algebra 0. Note that in this process, one can alternate the use of the 1-cocycle restrictions 
with the use of the O-coboundary cancellations. For instance, if the element p belongs to the 
algebra h, the function {F;p) = r/ p can be annihilated by the image of the function *$> p = j r] p dx 
under the O-coboundary map. Indeed (5o^ p ;p) = p(j r] p dx) = r\ p and F = F — 5o^ p belong 
to [F]. Thus, we can assume the function rj p to be equivalent to the zero function. Then for 
another vector field v in [), using the 1-cocycle restriction for the pair (p,v), that is 

p(F;v) -v(F;p) - (F;[p,v}) = prjv - - V\p,v] G 

one obtains conditions on the two functions r\ v and V[p,v}- Once again, one might try to absorb 
part of the function r\ v with 5o^> v , the image of another function fy v . Note that, in order to 
maintain rj p = 0, a restriction is imposed on fy v . Indeed, when the 1-cocycle F + 5q^> v is 
applied to p, we have to avoid reintroducing a function for rj p . Thus we need to consider only 
the functions *$> v for which (Sop; *$> v ) = (^ v ) x is a constant function. Then, to complete the 
determination of [F], the same process is preformed to every vector field of f), with some care in 
the choices of the O-coboundary maps, avoiding to undo the simplifications done in the previous 
steps. 

The results of this partial classification of cohomology classes [F], that gives a partial clas- 
sification of Lie algebras of differential operators g, are summarized in Tables 1 and 2 at the 
end of this section. The first table gives a 1-cocycle representative for the twelve type I Lie 
algebras and Table 2 exhibits the results for some general Lie algebras of vector fields among 
the type III Lie algebras. For these two tables, the classification numbers, given respectively 
by Lie and Amaldi, sit in the the first column. The second column gives a basis for the Lie 
algebra f) and the third column exhibits the first order differential operators v + rj v for which 
{F; v) = ij v in not trivial. rj v is taken to be the simplest representative and when the function 
r\ v is trivial, the differential operator is simply the vector field exhibited in the second column. 

It would be impractical to present the details of the computations in all cases. Furthermore, 
the arguments are quite similar for all Lie algebra f) of vector fields. So, for brevity's sake, we 
will only give the details for two of the selected Lie algebras. The chosen examples illustrate 
well the general process and will give to the reader a good idea of how the calculations proceed 
in general. 

Type I, case 1. This Lie algebra f) is spanned by the eight vector fields p, q, xp, yq, xq, 
yp, x 2 p + xyq and xyp + y 2 q. The vector field p belongs to the Lie algebra, hence, as mentioned 
previously, the function r] p can be assumed to be zero. The 1-cocycle condition for the pair (p, q) 
imposes the following restriction 

(5iF;p,q) = (rj q ) x - (r] p ) y - (F ; [p,q]) = (r) q ) x £ {1}. 

Thus rjg = c q x + h q (y, z), where c q is a constant. Hopefully, the function h q (y, z) can be absorbed 
by the image, under the O-coboundary map, of the function ^> q = j h q (y, z)dy. Since (^ q ) x is 
zero, the O-coboundary of ^ q will not affect 7] p . 
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Similarly, by considering the pair (p,xp), one concludes that r) xp = c xp x + h xp (y,z), where 
be canceled, without changing the previous functions, by the O-coboundary of the 
function ty xp = c xp x. Then, for the pair (q,xp), the restriction reads as 

(SiF; q, xp) = (r] xp ) y - x(r] q ) x - (F; [q, xp]) = (h xp (y, z)) y - x • c q E {1}. 

Necessarily, since h xp depends only on y and z, the constant c q has to be zero and the function 
h xp (y, z) is forced to be of the form d xp y + K(z), where d xp is a constant. Thus, at this point, 
Vp = °> Vg = and Vx P = d xp y + K(z). 

Consider now the three vector fields yq, xq and yp. If we pair each of them with p and q, 
from the six 1-cocycle restrictions, one obtains directly the following 

f]yq — CyqX -\- dy q y + k yq yZ), Vxq — c xq x + d xq y + k xq (z), T]y P — Cy p x + dy p y -\- kyp(z). 

With the image of the function ty yq = d yq y, the function rj yq can be reduced to rj yq = c yq x+k yq (z) 
without undoing the previous work. Prom the restriction associated to the pair (xp,yp), one 
easily check that 

{5 1 F; xp, yp) = x(r] yp ) x - y(j] xp ) x + r] yp = x- c yp + c yp x + d yp y + k yp (z) G {1}, 

forcing c yp and d yp to be zero and k yp (z) to be a constant function. Similarly, by considering 
the pair (xp,yq), one obtains that the function x ■ c yq — y ■ d xp must be constant, hence c yq 
and d xp are zero. To completely determine the functions r/ v for these three vector fields, two 
restrictions, associated to the pairs (xp,xq) and (yp,xq), must be verified. The first imposes 
that x(r] xq ) x — x(j] xp ) y — rj xq = x ■ c xq — c xq x — d xq y — k xq (z) must be constant. Thus it leaves no 
choice but to take d xq as the constant zero and k xq (z) as a constant function. Finally, the last 
restriction forces y(i] xq ) x — x{rj yp ) y — rj yq + ri xp = y ■ c xq — k yq (z) + K (z) to be a constant, hence c xq 
must be zero while k xq {z) must be equal, modulo the constant functions, to the function K(z). 
Putting together these restrictions, the image of the 1-cocycle F for the first six vector fields 
of f) can be described as r\ p = r\ q = r\ xq = r\ yp = and r\ xp = r\ yq = K(z). One easily checks that 
the remaining two restrictions are satisfied. 

To determine completely the 1-cocycle F, it remains to find its images for the two vector 
fields T := x 2 p + xyq and Q = xyp + y 2 q. For t]t, three restrictions are needed to reach that 
t]t = 3xK(z). Indeed, from the pair (p,T), the cocycle condition forces the following equality 
(vt)x — %Vxp — Vyq = C T, where ct is a constant. It is not hard to see that rjr must be equal 
to 3xK(z) + ctx + hx{y, z). From the pair (q, T), we get similarly that t\t = 3xK(z) + ctx + 
dry + kx(z). Finally, the restriction for the pair (xp, T) leads to 

x(rjr) x ~ T(rj xp ) - r] T = x ■ 3K(z) + x ■ c T - (SxK(z) + c T x + d T y + k T (z)) G {1}. 

Hence the constant dx dies out and kx{x) has to be a constant function. Note that, with these 
3 restrictions, rjr = 3x(K(z) + ct/3), but, by taking r\ xp = K(z) + cy/3, one gets the claimed 
result. By symmetry on x and y, the exact same arguments lead to t]q = 3yK(z). It is then 
straightforward to verify that the 1-cocycle F, given by the eight functions 

Vp = Vq = Vxq = Vyp = 0) Vx P = Vyq = K ( z ), T] P = 3xK(z) and r]Q = 3yK(z), 

satisfies all the other 1-cocycle conditions. Finally, the Lie algebra q associated to this triple 
(f), {1}, [F]) is the Lie algebra spanned by 

{p, q, xp + K(z),yp, xq, yq + K(z), x 2 p + xyq + 3xK(z),xyp + y 2 q + 3yK(z), 1} , 

where K(z) can be any function. 

Fortunately the calculations performed for a given Lie algebra f) can be repeated for any 
other Lie algebra sharing a subset of generators with fj. Note also that some ad hoc lemma's 
were used trough this work to simplify these calculations. For instance. 
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Lemma 1. Let i : M? — ► M 3 , (x,y) \—> (x,y,z) denote the inclusion map and suppose that 
l)o C r(i*TM 2 ), meaning that the generators of f)o depend on the variables x and y only. Let f) 
be a Lie algebra of vector fields on M 3 given by f) = f)o © {r, zr, z 2 r}. If, for non constant 
functions f(x,y) and g(x,y), the vector fields f(x,y)p and g(x,y)q belong to f) and if their 
associated images T)f(x,y)p an d Vg(x,y)q depend on x and y only, then 

H\t ) ,H^)/{l})=H\t)o,m 2 )/{l})®HH{r,zr,z 2 r},m)/{l})- 

Proof. Denote A := f(x, y)p and B := g{x, y)q. From the cocycle restrictions associated to the 
pairs (A, z l r), where % = 0, 1, 2, we obtain 

(6 1 F; A, zV) = A(r, zir ) - ^(r/ A ) 2 - (F; [A, z'r}) = f(x, y)( Vz * r ) x - - (F, 0) 

= f(x,y)(ri z i r ) x e {l}. 

Since f(x,y) is not constant, (rj z i r ) x must vanish, hence the functions n z i r depend on y and z. 
In a similar way, from the restrictions associated to the pairs (B,z' l r) it is straightforward to 
conclude that t] z % r = h l (z). Finally, for any element v in f)o, the function r\ v will depend on x 
and y only. Indeed, since r\ zr depends on z only, 

(S 1 F; v, zr) = v(n zr ) - z(j] v ) z - (F; [v, zr}) = - z(rj v ) z - (F, 0) = -z(rj v ) z G {1}. 

Therefore (r] v ) z must be zero, forcing the function r/ v to depend on x and y only. ■ 

Note that H 1 ({r,zr,z 2 r},J r (R)/{l}) is already well known. The 1-cocycle F associated the 
Lie algebra rj = {r, zr, z 2 r} is determined by three functions and the simplest representative is 
given by rj r = 0, n zr = and n z 2 r = dz, for any constant d. Thus, one can use this lemma to 
simplify some of the computations required in this classification problem. For instance, given f) 
the type I Lie algebra of vector fields given by case 10 in Table 1, the Lie algebra of differential 
operators g built from f) is obtained from a direct application of this lemma. 

Type I, case 10. The Lie algebra f) = Span{p, q, xp, yq, xq, yp, x 2 p+xyq, xyp+y 2 q, r, zr, z 2 r} 
can be decomposed as f)o © {r, zr, z 2 r} where f)o is the case 1 Lie algebra from the same table. It 
was shown in the previous calculations that the functions T] yp and n xq are zero, hence functions 
on x and y only. Thus the case 10 Lie algebra, along with its two vector fields yp and xq, 
satisfies the requirements of the Lemma [TJ Therefore, for the vector fields in the algebra f)o, 
the values of the 1-cocycle depend on x and y only, forcing K(z) to be c a constant function. It 
is then obvious that the 1-cocycle F is defined by eleven functions, were the three non-zero are 
given by t]t = cx, tjq = cy and r\ z i r = dx, for c and d any constants. The Lie algebra of first 
order differential operators g corresponding to this triple is then 

g = Span{p, q, xp, yq, xq, yp, x 2 p + xyq + cx, xyp + y 2 q + cy, r, zr, z 2 r + dz, 1}. 

It should be pointed here that H 1 ^, .F(M)/{1}) and H 1 (t),F(NL)) can also be determined 
alternatively using isomorphisms given in [20] and [21]. For the case that interests us, that is 
H 1 ^, ^(NVj/ll}), we fix a base point e and denote i the isotropy subalgebra. Provided the 
existence of a subalgebra a C rj which is complementary to i, one can show, see |20j . that 

iTH^(M)/{l})-# 2 0)/i). 

This isomorphism leads to an explicit method for constructing 1-cocycle representatives F in 
H 1 ^, JF(M)/{1}). We first choose a, a 2-cocycle representative of a class in H 2 (t)/i), and, for 
{v 1 , . . . , v n } a basis of rj, we denote a.ij = a(v l ,v^). If o = {v 1 , . . . , v m }, i = {v m+1 , . . . , v n }, 
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and c& are the structure constants of the Lie algebra f), a 1-cocycle in i? 1 (t),^ r (M)/{l}) will be 
obtained by solving first the following m(m — 1) equations 

(fj) ~ v j (fi)-J2 c %h = OLij , for 1 < i < j < m. 

k 

Once a non unique solution fx, . . . , f m is obtained, the remaining functions / m +i, . . . , f n are 
determined as the unique solution to the m(n — m) equations 

v% Uj) -v 3 {fi) -^Cijfk = aij, for 1 < i < m, m + 1 < j < n, 

k 

with initial conditions 

fi(e) = for m + 1 < i < n. 

Note however that this method can not be applied to all the three dimensional Lie algebras 
since the existence of the complementary Lie subalgebra is not guaranteed. For instance, the 
type III case 17A\ can not be treated using the isomorphism. Indeed, for the Lie algebra 

f) = {p, q, xp + zr, yq, x 2 p + (2x + az)zr, y 2 q}, 

and the base point e = (0,0,0), the isotropy algebra i is generated by the last four elements 
and the algebra a = {p, q} fails to be complementary, due to the absence of the element r in 
the Lie algebra. One can easily verify that Z 2 (f)/i) = {oi\ A 03, a\ A 05, 02 A 04, 02 A a$} and 
B 2 (t)/i) = {a.\ A 03,02 A 04}. One can observe at that point that the theorem does not hold, 
since the dimension of H 2 (l)/i) is two while the dimension of H 1 ^, ,F(M)/{1}) was computed to 
be three previously. Moreover applying the technique to the 2-cocycle a = c-a\ A 05 + d-a.2 AQ6, 
one gets a 1-cocycle that does not satisfies all the conditions that were not considered in the 
technique detailed above. 

2.3 Classification of quasi-exactly solvable Lie algebras 

of first order differential operators and the quantization condition 

The Lie algebras given in Tables 1 and 2 are the candidates for being quasi-exactly solvable 
Lie algebras, i.e. we might expect them to admit N a finite dimensional module of smooth 
functions J\f. In the investigation for these explicit finite dimensional modules, some new re- 
strictions are imposed on the 1-cocycles F. Indeed, as for the quasi-exactly solvable Lie algebras 
in lower dimensions, it comes out that a finite dimensional module exists only if the values of 
the functions rj v are taken in a certain discrete set. For this reason, this restriction is named 
quantization condition. The quasi-exactly solvable Lie algebras and their fixed modules can be 
found in Tables 3 and 4 for, respectively, the type I and the selected type III Lie algebras. The 
first column use the same classification numbers as in Tables 1 and 2 and a representative for 
the non-trivial quantized 1-cocycles is exhibited in the second column. Finally, N , the finite 
dimensional g-modules of functions are described in the last column. Once again the detailed 
calculations are repetitive and the essence of the work can be grasped with one or two examples, 
together with the following general principles. 

1. A finite dimensional module for the trivial Lie algebra g = {p} is defined as an x-translation 
module. For instance, any space spanned by a finite set of functions of the form 

n 

h = ^2g i {y,z)x i , 

i=0 
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along with all their x derivatives, is an x-translation module. This particular case of 
x-translation module is referred as a semi-polynomial x-translation module. The most 
general x-translation module is obtained by a direct sum 

AA = 0A/- A , N\ = N\e Xx , 
aga 

where M\ are semi-polynomial x-translation modules and the exponents are taken in 
a finite set A, read [9] for more details. Obviously, the y, and the ^-translation mod- 
ules are defined the exact same way. 

2. If the Lie algebra under consideration contains the two differential operators p and xp, 
the module M will be an x-translation module and the operator xp will impose extra con- 
straints. Firstly, all the exponents A need be zero. Otherwise, for an non-zero exponent A, 
the degree in x of the generating functions in the module M\ would be unbounded, con- 
tradicting the finite dimensionality of J\f. Moreover, if h = ^2™=$ g l (y, z)x l belongs to the 
module jV, the function xh x also needs to belong to that module. Note that both functions 
have the same degree in x and are linearly independent if h is not a monomial. Thus, by 
an appropriate linear combination of these two functions, one can reduce the number of 
summands in h. By iterating this process, each generating function can be reduced to 
a monomial in x. Thus, h = g(y, z)x l where g(y, z) belongs to G l a finite set of functions 
in y and z. Since M is a x-translation module, h x = ig(y, z)x % is also a function in A/", 
hence g(y,z) needs to be also contained in G' l ~ l . Therefore, the module J\f decomposes 
into the following direct sum 

M'=@x i gtiv,z), i = 0,...,n, k = 0,...,k, 

where all the functions gt(y, z) belong to G 1 a finite set and where G % C G i_1 . 

3. Likewise, if a Lie algebra g contains the elements p, q, xp, and yq, a general finite dimen- 
sional g-module for this Lie algebra will be at most 

N = (B xt y J 9k J ( z ), i = 0,...,n, j' = 0,...,m, k = 0, . . . ,l {iJ) , 
where the functions g 1 ^ (z) belong to G^\ a finite set of functions of z satisfying G^'^ C 

Gt'-ijlnGty'- 1 ). 

A simple method to describe these modules is to represent each generating function 
x l y 3 g 1 ' 3 (z) by a point (i, j), in the Cartesian plane. If a vertex belongs to the dia- 
gram, since M is an xy-translation module, the vertices (i — and (i,j — 1) must also 
sit in the diagram. To complete the description, a finite set G^'i' is associated to each of 
these vertices, with the same restriction as above. For instance, such module N can be 
represented by 

j 
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with all the sets G™' being equal to {z,e z }, with the exception of G^ 3,1 ) that contains 
only the function z. It is then straightforward to verify that this module is indeed 



, x 2 z, x 2 e z , xyz, xye 2 



y 2 z, 



N = Span{0, z, e z , xz, xe z , yz, ye z 

2 z 3 3 z 2 2z 2 1 z "\ 3 z 3 i 

y e ,x z,x e ,x yz,x ye ,xy z,xy e ,y z,y e ,x yz}, 



and that it is a fl-module for the Lie algebra q = {p, xp, y, yq}. 

If the Lie algebra q contains the differential operators p, q, xp, yq and yp, from the three 
previous principles, the generators for a g-module are given by h = x'y- 7 g 1 ^ (z). After 
applying the operator yp on h, the resulting function reads as ix yi + (z) . Thus, 
iterating this operator, we conclude that all the functions x l ~ r y :,+r g l,:> (z) must belong 
to N , for r < i. Since p and q also belong to the algebra, all the functions x a y h g % ^ {%), z) 
with a < i, b < j and a + b < c must belong to M. This condition can be expressed by 
the following inclusion G^ C G^~ 1 ^ n G^ j ~^ n G^~ l ^ +l \ and observe that the first 
set can be omitted without affecting the condition. To summarize, the g-module 

will be at most 



0, 



, n, 



0, 



,771, 



k 



0, ■ • • , 



where the functions g 1 ^ (z) belong to G^'^ , a finite set of functions with C n 

Once again it is possible to represent such module by a diagram along with a set of 
functions for each vertex of the diagram. The restrictions for these sets are G™' C 

q( 1 ~ 1 J+ 1 ) p| Q(i,j—1) anc j conditions on the vertices are slightly different from the one 
in the previous example. Indeed, if a vertex belongs to the diagram, the two vertices 

(i — l,j + 1) and — 1) must also belong to the diagram. Note again that this implies 
that the vertex (i — also lies in the diagram. For instance the diagram, 

j 



together with twenty appropriate sets of functions C?(w for each vertex, would generate a 
g-module for the algebra g = {p, q, xp, yq, yp}. 

5. Finally, if the elements p, q, xp, yq, yp and xq sit in the Lie algebra under consideration, 
the module will be at most 

AA = 0xV^'(z), i + j = 0,...,n, k = 0,...,l id , 

where the functions g l u 3 (z) belong to G^ %+ ^ , a finite set of functions with C G 1 -' -1 ). 
Indeed, consider h = g l ' J (z), a generator of bi-degree i + j = c. Since xq[h] and yp[h] 
must also lie in M, all the functions x a y h g % ^{z) with a + b = c will belong to J\f. Hence 
g {iJ) e and, reciprocally, g^ b) G G^ j \ Thus for all pairs (a, b) with a + b = c, 
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the finite sets G^ a ' b ^ are identical and it is therefore well defined to pose G^ a ' b ^ = G^ a+b \ 
Obviously, since M is a xy-translation module, the following inclusions hold C G^ l ~ l \ 

For these modules, the possible diagrams are more restricted and have necessarily the 
shape of a staircase. Also, instead of assigning one set of functions to each vertex, such 
a set is coupled to all the vertices having same total degree i + j- For instance, the module 
represented by the diagram 

j 



will be completely determined after fixing six sets of function in z. Note that this choice 
must respect the inclusion G®(z) C G^~ l \z), for i = 1, . . . , 5. 



Table 1. Cohomology for the type I Lie algebras of vector fields, m = {!}. 





Generators 


Cocycles 




1 


{p, q, xp, yq, xq, yp, x 2 p + xyq, xyp + y 2 q} 


xp 4 


- K(z), x 2 p + xyq - 
K(z),xyp + y 2 q-\ 


\-xK(z), 

- yK(z) 


2 


{xq, xp - yq, yp, 

Z\z)p,...,Z\z)p,Z\z)q,...,Z\z)q} 









3 


{xq, xp - yq, yp, xp + yq, 
Z 1 (z)p,...,Z l (z)p, Z\z)q,...,Z l (z)q} 


xp 4 


yq + K(z) 




4 


{p, Q, xp, yq, xq, yp, x 2 p + xyq, xyp + y 2 q, r} 


2 

x p 


+- xyq + cx, xyp + 


y 2 q + cy 


5 


{xq, xp — yq, yp, z k e x ' z p, z k e x ' z q, r}, 
k <ni, I = 0, . . . , b 









6 


{xq, xp - yq, yp, xp + yq, 

z k e XlZ p, z k e XlZ q, r}, k < n u I = 0, . . . ,b 


xp 4 


yq + cz 




7 


{P, q, xp, yq, xq, yp, x 2 p + xy, q, xyp + y 2 q, r, zr) 


2 

x p 


+- xyq + cx, xyp + 


y 2 q + cy 


8 


{xq, xp - yq, yp, p, zp, . . . , z l p, 
q,zq,..., z l q, r, zr + a(xp + yq)} 









9 


{xq, xp - yq, yp, xp + yq, p, zp, . . . , z l p, 
q,zq,..., z l q, r, zr}, k < n h I = 0, . . . ,b 









10 


{P, q, xp, yq, xq, yp, 

x 2 p + xyq, xyp + y 2 q, r, zr, z 2 r } 


2 

x p 
z 2 r - 


+- xyq + cx, xyp + 
f dz 


y 2 q + cy, 


11 


{xq, xp - yq, yp, p, zp, . . . , z l p, q, zq, . . . , z l q, 
r, zr + 7}(xp + yq), z 2 r + az(xp + yq)} 


z 2 r - 


f az(xp + yq) + cz 


12 


{xq, xp - yq, yp, xp + yq, p, zp, . . . , z l p, 
q,zq,..., z l q, r, zr, z 2 r + az(xp + yq)} 


z 2 r - 


f azixp + yq) + cz 



Note here that b, I and n; are positive integers, a, c, d, k and A/ are arbitrary constants and 
Z x {z), . . . , Z l (z) and K(z) functions of z. 
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Table 2. Cohomology for some type III Lie algebras of vector fields, m = {1}. 





(-IpTi prfltnrQ 


VjULj Llto 






4A 


{p, yq, q,xq,...,x t q + r,..., x l q + (*)x t_ V, . . . , 
x"q + (Jx r, xp — tzr, yq + zr}, < t < s 









4C 


{q, xq,..., x s q, p, yq, xp, x l y n ~ b r, zr}, 
(J < b < n, / < £q + so 









4D 


|g, xg, . . . , x s q, p, yq, xp, r, zr, z z r} 


2 i 

z r + cz 






5A* 


{q + r, xq + xr, . . . , x s g + x s r, p, 
xp, yq + zr, x 2 p + sxyq + sxzr} 


2 i 

x p + sxyq - 


- sr + cx 




5C 


{ S 2 I i / 7 I ~ \ 

{q,xq,. . . , x b q, p, yq, xp, x'p + sxyq + (/ + sn)xzr, 
x i y n-b r ^ zr ^ o < b < n, I < l + sb 


2 i 

x p + sxyq - 


1- (t + sn)r 


+ cx 


5D 


f s 2 i 2 1 

{q, xq, . . . , x q, p, yq, xp, x + sxyq, r, zr, z r\ 


2 i 

x p + sxyq - 


2 i 

- cx, z r + 


dz 


7C 


{p, 2xp + yq, x 2 p + xyq, x l y~ n r, zr}, < I < n 


x 2 p + xyq + 


cy 2 




17Ai 


{p, q, xp + zr, yq, x 2 p + (2x + az)zr, y 2 q} 


x 2 p + (2x + 
y q -f ay 


az)zr + bx 


+ cz, 


17A 2 


{p, q, xp + azr, yq + zr, x 2 p + laxzr, y 2 q + 2yzr} 


x 2 p + 2axzr 
y 2 q + 2yzr - 


+ cx, 
- dy 




17C 


{P, q, xp, yq, x 2 p + l xzr, y 2 q + p yzr, x l y p r}, 
l<l ,b<b 


x 2 p + loxzr 
y 2 q + b yzr 


4- cx, 
+ dy 




17D 


{p, xp, x 2 p, q, yp, y 2 q, r, zr, z 2 r} 


x 2 p + ax, y 2 


q + cy, z 2 r 


+ dz 



Note here that b, bo, I, lo, n, s and t are positive integers, a, c and d are arbitrary constants. Remark 
that Amaldi's Lie algebra 5A is not a Lie algebra. Indeed, for the space spanned by {q,xq, . . . ,x t q + 
r, . . . , x tJfl q + ( ()xV, . . . , X s q + ( s t )x s ~ t r,p, yp, xp — tzr, yq + zr, x 2 p + sxyq+ (s — 2i)xzr, x l y n ~ h r, zr} to 
be a Lie algebra, the parameter t needs to be zero. We then get the Lie algebra 5A* given in the table. 

This set of principles is of great help in the determination of the possible g-modules N for 
each Lie algebras of first order differential operators g described in Tables 1 and 2. Depending 
on the elements contained in the Lie algebra studied, we started our search of g-module based 
on the general module given in this guideline. Once again, the computations are tedious and it 
would not be relevant to detail each of them. We will concentrate on the same Lie algebras as 
in the previous step of this classification problem, that is the type I Lie algebras cases 1 and 10. 

Type I, case 1. Since the Lie algebra contains the differential operators p, q, xp, and yp, 
from the principle ([3]), the most general module M will be spanned by functions of the form h = 
x l yj g 1 '! (z), where g i '^{z) belongs to G^'^. We now consider the operator T = x 2 p+xyq+3xK(z) 
in the algebra q and its action on h = x n y a g n ' a (z) a generator of N with maximal exponent 
in x. Thus 

T[h] = nx n+1 y a g n ' a {z) + ax n+1 y a g n ' a (z) + 3K(z)x n+1 y a g n ' a (z) 
= (n + a + 3K(z))x n+1 y a g n ' a (z). 

Since the exponent in x was taken to be maximal, this imposes that K(z) is indeed a constant K 
equal to — Symmetrically, by considering Q = xyp + y 2 q + 3yK(z) and h = x b y m g b,m {z), 
a function with maximal exponent in y, the following equality holds K = — ^y^. For this to 
be possible, we necessarily have n + a = b + m. Consequently, the differential operators xq 
and yp belong to the Lie algebra g and the module J\f is given by the principle ([5]). Also note 
that the operators xp and yq force both a and b to be zero. Otherwise x n+1 y a ~ 1 g n ' a (z) and 
x b-iym+igb,m^ z ^ wou y [ n j\f ^ contradicting the maximality of n and m. Thus 3K(z) = —n 
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Table 3. Type I quasi-exactly solvable Lie algebras of differential operators and their fixed modules. 





Quantization condition 




Fixed module 




1 


2 i i 2 

x p + xyq — nx, xyp + y q — 


ny 


{x ■ y J g{z)\ 1+3 <n, g(z) <E G {l+J >\, 
for G« finite with G« C G ( ' _1) 




2 







{x*^g(z)\i + j<n, g(z)€G^ }, 
for G (l) finite with Z % G {1) C G ( ^ 1} 




3 







{x l y J g(z)\ i+j < n, g(z) G G^ + - 7 - ) }, 
for G w finite, and Z ? G m C Z^G^ 1 ) 




4 


x 2 p + xyq — nx, xyp + y 2 q — 


ny 


{x l y j z k e XlZ \ i + j < n, k < mi, 1 = 0,... 


M 


5 







{x l yiz k e x ' z \ i + j < n, z k e x ' z € G^}, 
for G y > a imite z-translation module, 
and Z l G {l) C G^ 1 ) for = z fc e A ' z p 




6 







{x l ?/-?,z fe e Ai2 i + j < n, k < mi, 1 = 0,... 
for G^ a finite ^-translation module, 
and Z l G^ C G^ 1 ) for Z l = z fc e Ai2 p 




7 


x 2 p + xyq — nx, xyp + y 2 q — 


ny 


{x l y j z k \ i + j < n,k < m} 




8 







{x^yi z k \ i + j < n, l(i + j) + k < m} 




9 







{x l y J z k \ i + j < n, l(i + j) + k < m} 




10 


x 2 p + xyq — nx, xyp + y 2 q — 
z 2 r — mz 


ny, 


{x l yiz k \ i + j < n,k < m} 




11 


z 2 r + az(xp + yq) — mz 




{x l y j z k \ i + j < n, a(i + j) + k < m} 




12 


z 2 r + az(xp + yq) — mz 




\x % y^z k \ i + j < n, a(i + j) + k < m} 





Note here that m and n are positive integers. 



Table 4. Some of the type III quasi-exactly solvable Lie algebras of differential operators and their fixed 
modules. 





Quantization condition 




Fixed module 


4A 







{x i y j z k 


i + sj + (s — t)k < n,j < m y , k < m z } 


4C 







{x l y J z k 
with bk- 


i + sj + (Iq + sn)k < n,k < b, j < b^} 
-l > h + n 


4D 


z 2 r — mz 




{x l y° z k 


i + sj < n,k < m} 


5A* 


x 2 p + sxyq + sr — nx 




{x l y j z k 


i + s(j + k) < n, j < m y , k < m z } 


5C 


x 2 p + sxyq + (Iq + sh)xzr — nx 




{x l y° z k 
with bk- 


i + sj + (Iq + sh)k < n,k < b, j < b^} 
-i>b k + h 


5D 


x 2 p + sxyq — nx, z 2 — mz 




{x l y 3 ' z k 


i + s j l£ n,k < m} 


7C 












17Ai 


x 2 p + (2x + az)zr — m x x, y 2 q — 


m y y 


{x i yi\ i 


<m x ,j< m y } 


17A 2 


x 2 p + 2axzr — m x x, y 2 q + 2yzr 


- m y y 


{x l y j z k 


i + 2ak < m x ,j + 2k < m y , k < m z } 


17C 


x 2 p + l Q xzr — m x x, y 2 q + b yzr 


- m v y 


{x l y° z k 


i + l a k < m x ,j + b k < m y , k < m z } 


17D 


x 2 p — m x x, y 2 q — m y y, z 2 r — m 


zZ 


{x % y J z k 


i < m x ,j < m y , k < m z } 



Note here that m, n, m x , m y and m z are positive integers. 
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and the module 

M = {./••//'//••'( •! \i + j<n, g l > j (z) e G^}, where C (8) 

is fixed by all the differential operators in g. Therefore it is possible conclude that the Lie 
algebra 

= Span{p, q, xp, yq, xq, yp, x 2 p + xyq - nx, xyp + y 2 q - ny, 1}, 

is quasi-exactly solvable with respect to the finite dimensional g-module N . 

Type I, case 10. Since the case 10 Lie algebra contains the case 1 Lie algebra, its module M 
will be at most the module given in ([5]). Observe first that the constant c in the case 10 Lie 
algebra has to be the negative integer —n. Furthermore, the operator r imposes N to be a z- 
translation module and the operator zr forces G^ to be generated by monomials. Then, for z m 
a monomial of maximal degree in , the function h = x l y l ~ % z m belongs to M . Since z 2 r + dz 
belongs to the Lie algebra, 

z 2 r + dz[h] = mxV- l z m+1 + dx^^z^ 1 = [m + d]x' V~' i z m+l , 

should belong to the g-module M. Thus, from the maximality of the degree in z, the constant d 
has to be the negative integer — m. Since the argument must hold for every set G^\ they will all 
share the same monomial of maximal degree m. We can therefore conclude that the Lie algebra 

g = Span{p, q, xp, yq, xq, yp, x 2 p + xyq — nx, xyp + y 2 q — ny, r, zr, z 2 r — mz, 1}, 

is quasi-exactly solvable with respect to the module 

M = {x l y j z k | i + j < n,k < m}. 

To summarize, a partial classification of quasi-exactly solvable Lie algebras of first order 
differential operators was accomplished in this section and the description of these Lie algebras g, 
along with their g-modules, can be found in Tables 1-4. In principle, it would be possible to 
achieve a complete classification using similar arguments. However this gigantic work would 
require a colossal amount of time. Nevertheless this partial classification is a good starting 
point for seeking new quasi-exactly solvable Schrodinger operators in three dimensions. In that 
scope, the next section is devoted to the description of few new quasi-exactly solvable Schrodinger 
operators. 

3 New quasi-exactly solvable Schrodinger operators 
in three dimensions 

Recall that in the general classification problem, once the quasi-exactly solvable Lie algebras of 
differential operators g are determined, the next step is to construct second order differential 
operators Ti that are locally equivalent to Schrodinger operators. Given g, one of the Lie algebras 
of first order differential operators obtained in the previous section, we obtain a second order 
differential operator Ti by letting 

m m 

H = C ab T a T b + ^2 C aT a + C , where T a G g, (9) 

a,b=l a=l 

as illustrated previously. Then, one has to choose the coefficients C a b, C a , Cq in such a way 
that the closure conditions du = are satisfied. Given T a = ^ ai di + n a , the closure conditions 
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are the Frobenius compatibility conditions for the overdetermined system 



a 



ah 



2rf 



C a 



0. 



where a = A + ^ ln(g) and fi = e A is the gauge factor. Finally, we need to bear in mind that the 
last step in the classification is to verify that the operators are normalizable, i.e. the functions 
in M, the module obtained after the gauge transformation, need to be square integrable. These 
operators will therefore have the property that part of their spectrum can be explicitly computed. 

We have now in hand a large variety of generating quasi-exactly solvable Lie algebras q. 
The door is therefore wide open to the construction of numerous new quasi-exactly solvable 
Schrodinger operators in three dimensions. However the Schrodinger operators described in this 
paper are built only from two of these new quasi-exactly solvable Lie algebras: the type III 
cases 17D and case 5A*. The reader can therefore see that many more examples can be 
constructed using this method together with the results of the previous section. 



3.1 Type III, case 17D, (sl(z) X sl(z) x sl(z)) 

The first two families of normalizable quasi-exactly solvable Schrodinger operators displayed in 
this section are similar to the operator given in the example (pQ). However, these two examples 
are more general. Indeed, for these two operators, the type III case 17D quasi-exactly solvable 
Lie algebra g is spanned by the first order differential operators 

p, xp, x 2 p - m x x, q, yq, y 2 q - m y y, r, zr, z 2 r - m z z, 

where m x , m y and m z are non negative integers and the module M mx m y m z is generated by the 
[m x + l)(rn y + l)(m z + 1) monomials 

x l y3 z k where < % < m x , j < m y and k < m z . 



3.1.1 First example 

The first family of operators is constructed with the following choice of coefficients 



( A 





1 

















1 \ 





B 

















1 





1 





c 











1 

















2A 





























2B 





























2C 

















1 











.4 





1 





1 

















B 





V i 

















1 





c J 



C c = [0, 0, 0, — 2Am x , — 2Bm y , — 2Cm z , 0, 0, 0] and Cq = —Am x — Bm y — Cm z . 

From the order two terms of these second order differential operators, the induced contravariant 
metric 

g(ij) i S 

obtained and reads as 

A{x 2 + l) 2 {x 2 + l){z 2 + 1) \ 

By A + 2{B + l)y 2 + B . (10) 

(x 2 + l)(z 2 + l) C(z 2 + 1) J 
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The determinant of the matrix is g = (1 - AC)(x 2 + I) 2 (By 4 + 2(B + l)y 2 + B)(z 2 + l) 2 and one 
easily verifies, for A, B and C positive and AB > 1, that the matrix is positive definite on R 3 . 
The operator can therefore be written as 

-2H = A + V + U, 

where A is the Laplace-Beltrami operator related to the metric (jlOp and where 

V = -2(x 2 + l){Am x x + m z z)p - 2m y (By 3 + By + y)q - 2(z 2 + \){Cm z z + m x x)r. 

From a direct computation, the closure conditions are verified and the gauge factor required to 
gauge transform 7i into a Schrodinger operator 7i$ is given by 

fj, = (x 2 + l)^(By 4 + 2(B + l)y 2 + B)=^ (z 2 + l)=r*. 
Once the transformation is performed, the equivalent operator reads as 
-2H = A + U, 

where the potential of the Schrodinger operator — % is given by a rational function in y. 

Note that the same three factors arise in both fi and g. This will simplify our computations 
while testing the square integrability of the functions in M. Indeed, a function in M is given by 
h = fix l y^ z k where the exponents i, j and k are smaller or equal to m x , m y and m z respectively. 
Our aim is to show that the triple integral 

JJM 2 



{jix l y 3 z ) ^fgdxdydz 

is finite. Obviously, it is sufficient to show the convergence of this integral for the monomials 
of maximal exponent. We can therefore focus on 

x 2m x y2m y z 2m z 

dxdydz. 



m v 1 

R3 (3.2 + i)nta+i(jgj^ + 2(B + l)y 2 + B) — + *(z 2 + 



Using Fubini's theorem, this triple integral can be factored into the product of three integrals 



x 2m x rco y 2m y 



dx, / my i dy, and 



(x 2 + l) m -+ 1 ' J-oc (By* + 2{B + l)y 2 +B)^^ 



oo z 2m z 
-oo (Z 2 + 1)^ +1 ^ 

each of which is easily shown to be convergent. We therefore have in hand a normalizable 
quasi-exactly solvable Schrodinger operator and it is feasible to determine explicitly part of its 
spectrum. 

For instance, if we fix m x = 0, m y = 2 and m z = 1, few manipulations lead to the six 
eigenfunctions of the operator 7i restricted to N . Indeed, with this choice of parameters, the 
g-module is 

= {i,y,y 2 ,z,yz,y 2 z}, 

and the transformation matrix to be diagonalized reads as 





-2 





2 


B 





o \ 







-A-2B 
















2B 





-2 






















-2 





2B 
















-A-2B 





V 











2B 





"2 / 
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Once the diagonalization is performed, three different eigenvalues Ai = —4 — 21?, A2 = —2 — 25, 
and A3 = —2 + 2B are obtained, each of them having multiplicity two. The six eigenfunctions 
are respectively 

ipi,i = y, ipi,2 = yz, ^2,i = i + y 2 , ^2,2 = -1 + y 2 , 
ip3,i = -z + y 2 z, ^3,2 = z + y 2 z. 

Consequently, we obtain three multiplicity two eigenvalues of the Schrodinger operator Tlo: 
Ai = 1 — B, X2 = 1 + B and A3 = 2 + B, and the six scaled eigenfunctions are 

ipi,i = fj>y, ipi,2 = fiyz, tp2,i = a*(i + y 2 ), ^2,2 = ^(-1 + y 2 ), 

^3,1 = V(~z + y 2 z), ^3,2 = n{z + y 2 z). 

As mentioned previously, the metric (|10p is positive definite on R 3 , hence Riemannian, and 
one can compute that the Riemann curvature tensor is zero everywhere. The change of variables 
that leads to a Cartesian coordinate system is given by 

f 1 

X = arctan x, Y= — , Z = arctan z, 

J v / By* + 2(B + l)y 2 + B 

where we have some flexibility on B to adjust the roots of the elliptic integral. Note that this 
operator is a good illustration of the modified Turbiner's conjecture in three dimensions proved 
in [6]. Indeed, the generating Lie algebra q is imprimitive and the leaves of the foliation are 
surfaces, the Riemann curvature tensor is zero and, as expected, the potential is separable since 
it depends on only one variable. 



3.1.2 Second example 

With the same representation of Lie algebra by first order differential operators but a different 
choice of coefficients, one constructs another family of second order differential operators H. 
Indeed, with 



1 A 





1 

















A \ 





D 























1 





B 











1 

















2A 





























(3D + C 





























2\B 

















1 











.4 





A 























0C 





V A 

















A 





X 2 B J 



C c = [0,0,0, -2Am x ,-0D(l + 2m y ) + C,-2XBm z , 0,0,0] and 
Co = —Am x — Bniy — Cm z , 

a family of operators 7i is obtained and one easily verifies that all these operators are equivalent 
to Schrodinger operators TCq. Note that this family of operators is slightly more general than the 
family obtained in the first example. However, some of the details are lengthy and are omitted 
for brevity sake. The induced contravariant metric is given by: 

A{x 2 + l) 2 {x 2 + l){\z 2 + 1) \ 

0Cy A + y 2 {2(5 + 0D + C) + D , (11) 

(x 2 + l){Xz 2 + 1) B{Xz 2 + z)) ) 



Quasi-Exactly Solvable Schrodinger Operators in Three Dimensions 



21 



and it is positive definite on M 3 provided A, B, C, D and (5 are positive and AB > 1. Its 
determinant is g = (AB - 1) (x 2 + l) 2 (/3Cy 4 + 2f3y 2 + (3Dy 2 + y 2 C + D) ( Xz 2 + l) 2 , and the gauge 
factor required is the product of three functions in x, y and z respectively. After the gauge 
transformation the new potential is again a rational function involving only the y variable and 
the Riemann curvature tensor is null again. The following change of variables leads to Cartesian 
coordinates 



X = arctan x, Y = J — ; 



Z 



1 

+ 2/V + 2f3Dy 2 + y 2 c + D 

arctan \f\~z 



where we have some flexibility on f3, C and D to adjust the roots of the elliptic integral Y. 

However, we do not know if these operators are all normalizable. But, if we fix C = (3D, the 
gauge transformation simplifies and becomes, once again, very similar to the determinant of the 
metric (1111). Indeed 



^ = (x 2 + 1)-T(/3 2 IV + 2/% 2 (l + D) + DY^{\z 2 + 1)"T 

and one verifies, the exact same way as in the previous example, that the functions in are 
square integrable. Therefore, for any choice of integers m x , m y , and m z , one would obtain 
(m x + l)(m y + l)(m z + 1) eigenfunctions in the spectrum of the Schrodinger operator TCq. 

For instance, if we fix A = 1, (3 = 5 and the three parameters m x , m y and m z to be 1, one 
gets two eigenvalues, —3 and —7 of multiplicity four, and the following eight eigenfunctions 

^-7,1 =^(-1 + xz), Tp- 7t2 = n(y - xyz), 

tp-7,3 = v(xy + yz), ^-7,4 = n(x + z), 
= (J>(y + xyz), ^-3,2 = K~ x + z )i 
= v(-xy + yz), ip-3,4 = m(i + xz). 

Note that the nodal surfaces can described easily in this coordinate system. Indeed, since (i 
is always positive, the nodal surfaces are simply the zero loci of polynomials. For these eight 
eigenfunctions, the surfaces are given by the zeros of degree two factorizable polynomials and 
one easily gets the following pictures. 
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xy-yz 

xy+yz 




3.2 Type III, case 5A*, sl(z) x C s+1 

For the last example, we consider the type III case 5A* quasi-exactly solvable Lie algebra and 
we fix the parameter s to be one. This Lie algebra q is therefore spanned by the following six 
first order differential operators 

p, q + r,xq + xr, xp, yq + zr, and x 2 p + xyq + xzr — nx, 

and from the Table 4, the g-module of function is given by 

M = {x l yi z k \ i + j + k < n,j < m y , k < m z }, 
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where n, m y and m z are non-negative integers. A family of Schrodinger operators on R 3 \{x = y} 
is obtained from the following choice of coefficients, 



/^0 0000\ 

B 

C 



D 

\ 0000/ 

C c = [0,0,(7,0,-2(1 + m)D,0], and 



ab 



C = (l + m) 2 D, 



where the parameters A, B, C, and D are positive. The induced contravariant metric is given 
by 



Cx 2 + A 
Dy 2 + B Dyz + B 
Dyz + B Dz 2 + B 



(12) 



its determinant is g = BD(Cx 2 + A)(y — z) 2 and the metric is positive definite on ]R 3 \{x = y}. 
Before performing the gauge transformation the operator reads as 

-2H = A + (2Cx - Cmx)p + (-Dy - 2Dmy)q + (-Dz - 2Dmz)r 
- l/2Cm + 1/ACm 2 + (1 + m) 2 D, 

and one easily verifies that the operator respects the closure condition. The gauge factor required 
to obtain a Schrodinger operator is 



H = (Cx 2 + A)—(y-z) 2 , 

and once again, contains the same factors as the determinant of the covariant metric. Finally, 
after the gauge transformation, the Schrodinger operator reads as, 

-2W = A + U, 

where U depends on the three variables. Although, it is not known if the functions in J\f are 
square integrable on the domain M 3 \{:r = y}. 

Note that for this example, the scalar curvature is constant and depends on the parameter D 
while the Riemann curvature tensor is equal to 



; dydzdydz. 



B(y - zf 

However, the potential does not seem to be separable. 
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